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We investigate the properties of electron states localized at the edge of a superconducting island 
placed on the surface of a topological insulator in a magnetic field. In such systems, Majorana 
fermions emerge if an odd number of vortices (or odd multivortex vorticity) is hosted by the island; 
otherwise, no Majorana states exist. Majorana states emerge in pairs: one state is localized near the 
vortex core, and another at the island edge. We analyze in detail the robustness of Majorana fermions 
at the edge of the island threaded by a single vortex. If the system parameters are optimized, the 
energy gap between the Majorana fermion and the first excited state at the edge is of the order of 
the superconducting gap induced on the surface of the topological insulator. The stability of the 
Majorana fermion state against a variation of the gate voltage and its sensitivity to the magnetic 
field allows one to experimentally distinguish the edge Majorana fermion from conventional Dirac 
fermions. 

PACS numbers: 71.10. Pm, 03.67. Lx, 74.45.+C 


I. INTRODUCTION 

The possible realization of Majorana fermion states in 
condensed matter physics is attracting considerable in¬ 
terest in recent yearsic—. This is partly due to the non- 
abelian anyonic statistics of Majorana fermions, allow¬ 
ing the realization of topologically-protected quantum 
gates^. Topological quantum computation requires the 
braiding of anyons^^. Majorana braiding might be real¬ 
ized by the controllable manipulation of the pairwise in¬ 
teraction between separate Majorana fermionsi^c—. The 
decoherence caused by the tunneling between Majorana 
fermions sets an upper limit of the time on the elementary 
operation, while the energy of the excited states deter¬ 
mines the lower time limit of the elementary operationi^. 
Many attempts have been performed to find Majorana 
fermions in different systems. Recently, possible obser¬ 
vations of Majorana fermions in quantum nanowire s^^d^ 
and atom chainsii were reported. 

The interface between a topological insulator and an 
s-wave superconductor (SC) is a promising system for 
the possible realization of Majorana fermionsJ^^— Such 
proximity-induced superconductivity is a mixture of s 
and p-wave correlationsBeing topologically equiv¬ 
alent to the p-wave superconductivity, it supports Ma¬ 
jorana fermionsj^ However, specific features of this su¬ 
perconducting state, in particular, quasiparticle linear 
dispersion on the surface of the topological insulator, 
are of importance for the structure and robustness of 
the Majorana state and requires adequate analysis. Ma¬ 
jorana fermions may emerge at the vortex core in the 
proximity-induced superconducting region on the sur¬ 
face of the topological insulator. However, the mini¬ 
gap separating the Majorana fermion and the Caroli-de 


Gennes-Matricon (CdGM) levels in the Abrikosov vor¬ 
tex cor o^^d^ is very small (about 10“^ K). In order to 
increase the robustness of the Majorana state. Ref. [2l| 
proposed to make a hole in the superconducting layer to 
pin the vortex and to remove the CdGM levels. These 
ideas were further elaborated in Ref. [13. 

Majorana fermions could also localize near boundaries 
of the superconducting region. However, such states 
have attracted much less attention than the Majorana 
fermions near the vortex core. In Refs. 131^ it has been 
shown that the edge Majorana fermion localizes at the 
interface between p-wave superconductor and a topologi¬ 
cally trivial insulator if odd number of vortices penetrates 
superconductor. Edge Majorana modes were studied the- 
oretically^S in a finite-size heterostructure made of a SC, 
a ferromagnetic insulator, and semiconductor with strong 
spin-orbit coupling. It has been argued in Refs. Mm 
that Majorana fermions can arise at the edge of a semi- 
finite SC placed on the surface of a topological insulator 
in a magnetic field. 

Significant progress with making superconducting is¬ 
lands on the surface of insulators or metals has been 
achieved in recent years. Observations of vortices and 
multivortices in the Pb superconducting islands^^— and 
regular structures of the Nb superconducting islands^ 
have been reported. These systems are of special inter¬ 
est for possible implementation of the Majorana fermion 
surface codes for topological quantum computation s"^^ . 

In our previous works^i^^l we investigated the Majo¬ 
rana fermion in the core of the vortex pinned by a hollow 
channel in s-wave superconductor placed on the top of 
topological insulator. This channel removes GdGM lev¬ 
els in the core of the vortex in s-wave superconductor 
making Majorana fermion robust^. 

In this paper we consider different system: thin cylin- 
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where i/Ti describes the topological insulator (TI), Hsc 
describes the s-wave superconductor (SC). The term T 
accounts for the tunneling from the TI to the SC, while 
rl represents the opposite processes: tunneling from the 
SC to the TI. The corresponding Bogolyubov-de Gennes 
equations are (after setting fi = e = c = 1) 

HtiiPti + T'l'^sc = wV’Ti, (2) 


FIG. 1: (Color online) Proposed experimental setup for the 
detection of edge Majorana fermions. A superconducting 
(SC) island of radius R (blue) is placed on top of a topological 
insulator (yellow). An external magnetic field is perpendic¬ 
ular to the interface. The black circle at the center of the 
superconducting island represents a vortex. 


drical s-wave superconducting island of radius R placed 
on the infinite surface of a topological insulator in a trans¬ 
verse magnetic held B (see Fig. [T|). If the island traps a 
vortex, two Majorana fermions are induced, one at the 
vortex core in the center of the island and the other at 
the island edge. 

This paper studies the Majorana states localized near 
the edge of a superconducting island. Edge Majorana 
states exist only if the vorticity I threading the island is 
odd, and disappear if I is even or zero. The energy split¬ 
ting between edge and vortex core Majorana fermions 
decays exponentially, when increasing the radius of the 
island. We will demonstrate that the edge Majorana 
fermion is robust: the gap between the Majorana fermion 
state and the edge excited state is of the order of the in¬ 
duced superconducting gap. This makes the edge Majo¬ 
rana state promising for experimental observation, and, 
possibly, manipulation. 

The presentation below is organized as follows. In 
Sec. U we derive the Bogolyubov-de Gennes equations 
for our system. In Sec. IIIII we investigate modes with 
zero energy for different number of vortices in the island. 
In Sec. m the system with a single vortex is considered. 
In Sec. |V] the obtained results are discussed, and conclu¬ 
sions are presented. 


Hsc'fpsc + TV'ti = w^/>sc- (3) 

The terms Hti, and Hsc can be written as 4 x 4 matrices 
in the Nambu basis 


Hti 

Hsc 


iv{a ■ Vr) - U{r)]Tz + v{a ■ A)to, 
■ (VR-fAr^)2' _ 

- "Z - Tz 

2m 




( 4 ) 


and T = r). In these equations, R = {x,y,z) 

is a point in the bulk of the SG, r = (a;,y) is a point 
on the surface of the TI, aj , Tj are the spin and charge 
Pauli matrices, A', A" are the real and imaginary parts 
of the order parameter in the s-wave superconductor, v 
is the Fermi velocity of the electrons on the surface of 
the TI, E-p is the Fermi energy in the SG, U (r) is a gate 
voltage applied to control the Fermi level in the TI, and 
A is the vector potential of the magnetic field. The wave 
functions ■0 ti, sc are four-component spinors 


V’TI.SC = 


( 5 ) 


It is easy to check that H satisfies the following charge- 
conjugation symmetry condition 


{H,E} = 0, 

E = CTyTyK, ( 6 ) 


where K is the complex conjugation operator. Conse¬ 
quently, for every eigenstate tp of the Hamiltonian H 
with a non-zero eigenenergy w 0, an eigenstate Eip 
with eigenenergy —w exists. The latter symmetry is ro¬ 
bust: small disorder does not destroy this property. 


II. BOGOLYUBOV-DE GENNES EQUATIONS 
A. Microscopic model 

Our system is shown in Fig. [TJ It consists of a thin 
cylindrical superconducting island of radius R placed on 
the surface of a topological insulator. The entire het¬ 
erostructure is in a magnetic field perpendicular to the 
topological insulator surface. 

To study the microscopic properties of such a system 
we will use the formalism of Ref. [l^. The Hamiltonian 
can be written as 

( 1 ) 


B. Effective Hamiltonian 

In this subsection we will derive the effective descrip¬ 
tion for the wave function of the electrons on the surface 
of the TI. To this end we exclude tpsc from Eqs. © 
and ([3|) to obtain 

{Hti + SjV'Ti = wV'Ti, (7) 

S = rt(w-ifsc)-^T. (8) 

We are interested in the bound states with energies lying 
within the superconducting energy gap, |w| < |A|. In this 
case, the self-energy matrix E was calculated in Refs, [l^. 


H = Hti + Hsc + T + T^, 
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[2l|. For the low-lying electron states with wave vectors k 
near the Dirac cone apex M in the Brillouin zone of the 
TI, k M, the self-energy is equal to 


S 


]VI,a; 


^ At^; - UJTQ 

v'lAP-cc^ 


- SUtz 


(9) 


C. Equations for the effective wave function 

Further, we assume that the island radius R is much 
larger than the SC coherence length ^sc- 

R > ^sc (17) 


where tq is the 2x2 unit matrix. This expression for 
S is independent of A. Such an approximation is valid 
for weak magnetic fields, which is assumed to be the case 
everywhere in this paper. The parameter A has the di¬ 
mension of energy. It characterizes the transparency of 
the interface between the TI and the SC: when A ^ Ey 
(A <C £’f), the barrier is transparent (non-transparent). 
The value 5U = 0(A) is the shift of the chemical poten¬ 
tial of the TI due to doping by the carriers coming from 
the SC. 

Using Eq. ((91) we can rewrite the Bogolyubov-de Gen- 
nes equation (l7|) in the form 


We are looking for solutions of Eq. cni) which correspond 
to bound states. The energies u of such eigenstates are 
smaller than the proximity-induced gap Axi, in Eq. (I16|l . 
If an Abrikosov vortex with vorticity I is trapped in the 
island, the order parameter A(r) can be expressed as 

A(r) = |A(r)| exp (—1/0), (18) 

where r and 6 are the polar coordinates, and |A(r)| —>■ 
|A|, when r ^ If the island radius R is large, i? ^ 
|A(r)| can be approximated as 

|A(r)| = |A|0(i?-r), (19) 


HeSlpTl = wV'TI, (10) 

where the effective Hamiltonian HeS equals to 

Hes = [iv{u;){a-Vr) - U{uj)]tz+ v{u;){a ■ A)to 

+ A'{uj)Tx + h"{w)Ty. (II) 

The renormalized parameters of the effective Hamilto¬ 
nian are 


where 0(r) is the Heaviside step function. In the geom¬ 
etry shown in Fig. (TJ the vector potential can be written 
as Az = Ar = 0, Aq = A{r). This choice of the vector 
potential corresponds to the magnetic field 


I d{rA) 
r dr 


( 20 ) 


Let us introduce a spinor F 


v{uj) 


V|AP-cc2 + A’ 


(t/-f (5t/)\/|A|2-a;2 




v/|A|2-u;2 + A 


( 12 ) 


tpTi = exp[-i9{lTz-az)/2 + i^j,9]F>^{r), 

= (/l^/2^/3^-/4T■ (21) 


(13) 


Here /i is the total angular momentum of an eigenstate. 
The transformation Eq. (I21() is well-defined only when 

j = M ”1 7^ (22) 


A(a;) 


AA 

\/|A|2-a;2 + A' 


(14) 


It is also convenient to define the renormalized coherence 
length: 




uy^|A|2 — a;2 ^ i/jAp — a;2 

AA A 


(15) 


If the Hamiltonian parameters are independent of r, 
then the eigenstates of the effective Hamiltonian obey the 
inequality |a;| > Axi, where the quantity Axi satisfies the 
following equation: 


Axi /a — Axi 
A V A -|- Axi 


(16) 


The physical meaning of Axi is the proximity-induced 
superconducting gap on the surface of the TI. 


is an integer. In other words, when the number of vortices 
I is odd (even), the angular momentum /r is an integer 
(half-integer). 

Substituting Eqs. (fTTl) . (fTTl) . and (1211) in Eq. (nni) we 
derive 

—^+^(^)^/f-|A|/i‘-(w+C)/2^=0, 

+ V A(r))/|‘ + |A|/f-(u;-U)/3^=0, (23) 

—^-^(O^/3^-|A|/2^-(a;-U)/{‘=0. 


Equations (IMI) have the following symmetries: (i) /r e->- 
/4 O i/i, /3 O 1 / 2 , U O -U, and {ii) A{r) ^ 
-Mr), fi ^ / 2 , fs ^ -fi, Z /i O- Therefore, 

we can consider further only A{r) > 0. 
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III. ZERO-ENERGY SOLUTIONS 


(l24ll and (1251) becomes 


If w, [/ = 0, the system of Eqs. (1^51) decouples into two 
sets of equations 


*^(|: + ^^^^^+A(r))/r + |A|/r=0, (24) 

+ l^lfs = 0 , 


The parameter |A| is zero outside the SC island area ac¬ 
cording to Eqs. (fT^ and (IT4)l . It is also zero at the center 
of the vortex core. In this paper we are mainly interested 
in the states localized at the edge of the island and the 
details related to the states in the vortex core are not of 
importance here. Then, for simplicity, we approximate 
the vortex core of the vortex with vorticity I by a cylin¬ 
drical hole with a radius about the coherence length 


fi = Ai exp 


fi = Ai exp 


/2 = A 2 exp 


/a = As exp 


J A{r')dr' 
0 







(27) 


The functions /2 and /a diverge when r —>• -|-oo; then, 
^2 = Aa = 0 or /2 = /a = 0. Matching solutions at 
r = R, one completes the derivation of the wave function. 

These eigenfunctions correspond to the a; = 0 Landau 
level, whose states are weakly corrected to account for 
the presence of the superconducting island. They are 
not Majorana fermion states, and a weak perturbation 
of the Hamiltonian may shift their eigenenergies away 
from zero value. 


A. System without vortices 


B. System with vortices 


Let us assume first that there are no vortices in the 
island. The magnetic field localizes a zero mode near the 
edge of the island. However, we show that this mode is 
not a robust Majorana fermion. 

In the absence of vortices, the SC order parameter A is 
non-zero in the SC island and the solutions of Eqs. dm) 
and (1^ which are regular at r = 0 can be expressed in 
terms of the modified Bessel functions Im{x) 


/i =(71 exp I - / A{r')dr' |4-i/2|^—j , 

fi = -iCiexp ^ J A(r')dr'j4+i/2^^^ , 

/2 = (72exp j A(r')dr'j4+i/2^^^ , 

/3 = *C'2exp I J A(r')dr'j4_i/2^^^ . (26) 


In this subsection we study a system with vortices. 
Since we are mainly interested in edge states and a rel¬ 
atively small SC island, we assume that the magnetic 
flux captured in the SC island forms a multivortex with 
vorticity I and further assume that the SC order param¬ 
eter in this vortex behaves like a Heaviside step func¬ 
tion A(r) = A0(r — ^). This simplification neglects 
the CdCM states inside the superconductor, and thus, it 
can significantly affect the Majorana state localized near 
the vortex core. However, the edge-localized Majorana 
fermion is fairly insensitive to the details at the center of 
the island, because its wave function decays quickly away 
from the edge. 

Under these assumptions, the solutions of Eqs. (IMl) 
and (1^ outside of the island r > R, and inside of the 
vortex core r < ^, can be expressed as 


/i = Ciexp \ - J A{r')dr' | r'"+ 


- / A{r')dr' r— 


Outside the island, where A = 0, the solution of Eqs. 


fi = Ci exp 


(28) 
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and 


/2 = C 2 exp J A{r')dr'j r ^ ‘ 2 " ^ 

/3 = C3exp (jA{r')dr^ 

where Ci are constants (for r < ^ these constants may be 
different from the constants at r > i?). The functions /2 
and fs diverge when r —> 00 , and also /2 diverges when 
r — >■ 0. Then, as it follows from Eqs. (ESI), a regular 
solution exists only if /2 = /s = 0 in the whole space. 

Inside the SC island a wave function can be expressed 
as a sum of two distinctly dissimilar solutions of Eq. (l24ll . 
A solution of the first type is localized near the vortex 
core: 

/i = (72 exp ^ j A(r')dr'j , 

f 4 = iC 2 exp ^ J A(r')dr'j . (29) 

Here (7( 2 are constants, and Kn{x) is the modified Bessel 
function. Since diverges at a; = 0, the function can 
be normalized only if /i < (Z + l)/2. Using the symmetry 
between positive and negative p,, one can generalize this 
inequality for arbitrary p,: 


\d\ < 


1 + 1 


(30) 

does not define 


If /i violates this condition, then Eq. 
a valid eigenstate. 

Unlike Eq. EH), which describes eigenfunctions local¬ 
ized at r = 0, a solution of the second type grows toward 
the edge of the island: for r < i? one can write 

fi = C[exp ^ j A{r')dr^ > 

f 4 = -iC[exp ^ j A(r')dr'j , (31) 

where /„ is the modified Bessel function of the second 
kind. Outside the island (r > R) the eigenfunction is 
defined by Eq. (1^51) . Most of the wave function weight 
is localized away from the island center. The value of r, 
where the weight is concentrated, grows as |^| increases. 
In the limit |p,| —>■ 00 , the wave function is virtually un¬ 
affected by the presence of a superconducting island at 
the origin. Thus, for large |/r| the eigenstates described 
by Eq. (pil) become indistinguishable from the states be¬ 
longing to the w = 0 Landau level of the Dirac-Weyl 
electrons. 


Majorana fermion states correspond to ^ = 0 solutions. 
In finite systems these states appear in pairs. In our 
case, one Majorana fermion is localized near the origin 
and another is localized at the edge of the island. Inside 
the island their wave functions are given by Eq. (l29l) and 
Eq. (inj). Outside the island, Eq. (1^ must be used. 

To demonstrate the “Majorana nature” of the p, = 0 
solutions, let us calculate the first-order corrections to 
the energies of the states Eq. ED) caused by a non-zero, 
but small, U: 


G= / drr'" Re ° 

r>R 


Sloi^ = UCif,2^G, 

— 2 J A{r')dr' 


(32) 


■ 1 2/.t 


iRi LiiJ 


ri “ 2 /^ 


— 2 f A(r')dr' 

drrRe ° 


( 

r 

2/7 

r 


i. 






-2ti\ 


V^V2(^ 


(Al) r2 (+ry 
\ V J V J 


i<r<R 


(\R\ t2 /'Aflj / ’ 

-I/2V V ) ^ii+l/2\ V )J 


where Ci^ are wave function normalization constants. 
For states localized near the vortex core of the vortex, 
we have a similar result with the modihed Bessel func¬ 
tions Ifj,±i /2 replaced by the modified Bessel functions 
K^±i/2- 

It follows from Eq. (IMl) that, if /i = 0, the correc¬ 
tion vanishes identically. Moreover, one can demonstrate 
that the fi = 0 states are invariant under the action of 
the charge-conjugation operator S, Eq. (I^. Thus, these 
eigenstates are topologically-protected mixtures of elec¬ 
tron and hole states. 

For all other values of /r this correction is nonzero. 
Thus, the zero-energy states with fi ^ 0 are not 
topologically-protected from the local perturbations of 
the chemical potential. In particular, the robustness of 
the zero-energy modes against variations of the gate volt¬ 
age may be used to distinguish Majorana states from 
conventional Dirac fermions. 

Finally, we would like to remind that, since /r is an 
integer only for odd I, see Eq. (l22l) . we must have an odd 
number of vortices on the island to generate Majorana 
fermion states. 


IV. SINGLE VORTEX 


In this section we study a system with a single vor¬ 
tex, / = 1, which is the simplest in terms of a possible 
experimental realization. We start with the case (7 = 0. 
According to Eqs. (I5T]) the wave function of the zero- 
energy state with /i = 0 localized near the vortex core 
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is 


■ 7 r ~ f ^(r')dr' — f ^ dr' 

= Bye~^^e 0 e 0 ” 


i 


0 

0 ’ 
-1 


This condition suggests that the growth of the tunneling 
parameter A improves the isolation of the two Majorana 
fermions from each other, which is a desirable property 
(33) for reliable Majorana state detection. 


where we used the explicit expression for the modified 
Bessel function of half-integer order. For the edge zero- 
energy eigenstate, the wave function is [see Eqs. (ESI)] 

11 


—i_ 

Here By and Bg are normalizing coefficients, which we 
choose real. Applying the particle-hole-conjugation op¬ 
erator S, we obtain that 'Eijjy(^f,-^ = ipv(^e) by direct calcu¬ 
lations. Hence, ipy and ^pe are Majorana fermion wave 
functions. 


- TT ~ f A{r')dr' f ^ dr' 

ipe = e 0 e“ ” 


A. Zero-mode splitting 

The hybridization between the vortex core and the 
edge Majorana fermions can, in general, split the degen¬ 
erate zero level. The splitting is zero in the case U = 0 
(see also Ref. [ 13 ). If t/ 0, the degeneracy is lifted, and 
Majorana fermions at the vortex core and at the edge 
form the two usual Dirac states. The wave functions 
of these states can be written as 4’± = T *'0e)/2. 
These functions satisfy the particle-hole symmetry of the 
Bogolyubov-de Gennes equations, = '0_. If we de¬ 
note the splitting energy as 2E +, then, the wave function 
corresponds to A+, while ip- corresponds to —E+. 
Let us assume that U is small. The first-order contribu¬ 
tion to the energy splitting becomes 


B. Excited states 


Now we calculate the energies of the excited states 
localized near the SC island edge, assuming that both 
Eq. (15^1) and Eq. (1551) are valid, and 17 = 0. It is implied 
for simplicity that the magnetic field penetrates the is¬ 
land uniformly, so, A{r) = r/2l‘^. In this paper, the 
states localized near the vortex core^ii^SiSl are not dis¬ 
cussed, because, when the inequality Eq. (1551) is valid, 
the vortex core states do not mix with the edge states, 
and may be neglected. 

When w 0, 17 = 0, ^ = 1, and A = 0, the system 
Eq. (1251) decouples. As a result, for r > i? we have two 
independent systems of equations: one for the electron 
components /i 2, another for the hole component /a 4: 




dr 


/r- I 


2q 




(39) 


Note that after the transformation if^ ^ and 

fs if 2 ^ 1 tbe first two equations and the second two 
equations exchange places. Thus, we can solve only the 
equations for /i_ 2 - Substituting / 2 , one derives for fi: 


E+=U{iP+\t.\iP+). (35) 

It is reasonable to assume that the applied magnetic 
field is smaller than the upper critical field, that is, 
h ^ where lb is the magnetic length, which satisfies 
lb = in the units used here. In this case we can 

neglect the effect of the magnetic field on the wave func¬ 
tions near the vortex corei^ Further, if the SC island is 
not large: 

i?<4, (36) 

then, using Eqs. (155|) . (154)) and (1551) we derive an estimate 
for the energy shift E^ in the form 


77+ oc — exp 

m 


-77/1(0) 


(37) 


dr‘2 r dr ^ V 


^ + ^1=0. (40) 


Solutions of the latter equation can be expressed in terms 
of the Tricomi confluent hypergeometric functions, tra¬ 
ditionally denoted^ as t/(a, b, z) (do not confuse it with 
the shift of the chemical potential U = U{r), which is a 
function of a single variable). As a result we have 


fi = exp U 


ti)‘^l? r^ 


2{)2 


.(41) 


If r llujfv, these wave functions decay as follows: 


One can now see that Majorana fermion states are robust 
against chemical potential variations if the radius of the 
SC island is large, in the sense that 

7?»|(0) = ^. (38) 


/!- = .C exp (-Ej . /? = C /f exp ; (42) 

consequently, they are normalizable. The second linear- 
independent solution to Eq. (001) diverges when r —>■ 00 , 
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thus, it is not included. Near the edge (r Ri R) we can 
approximate the functions in Eq. 65) as 

/r = tCJ^ (^) , = CJ^+i , (43) 

if the condition jv ^ 1 is satisfied. Later, we will show 
that this condition is similar to the initial assumption 
lb 3> R. The asymptotic behavior given by Eq. 65 may 
be guessed from Eq. 65- Indeed, near the island edge, 
r > i?, the terms and {fj,+l)/ll are much smaller 

than the remaining two, and may be omitted. After this 
simplification, equation (HOI) transforms into the Bessel 
equation. 

In the region r < R, we can neglect tl^e effect of the 
vector potential since R<^ lb- We introduce the following 
linear combinations^i 



^f = */r+/4, = 

Yr=if^ + fS, 

where Xi ^2 obey the differential equations: 

_ f 1 .A.ll]x^ = n 

dr'i 7. dr y[C('^)]^ vr J ^ ’ 

IdX^ ( I A 

dr2 ^ r dr yC(w)]2 hr r2 j ^ “ 0’ 

V 

where ^(w) = — — 


(44) 


(45) 

,(46) 


^|A|2-a.2 

and Fi ,2 can be found according to the relations: 






iv 1 

' dX^ 

A „ 

= — 


-Af 

UJ \ 

y dr 

V ^ 

iv 

i dX!/ 


to 

( dr 

+ — 

V 



(47) 


As it follows from Eqs. (l23l) . which are regular at r = 0, 
the solutions for Ai _2 can be expressed in terms of the 
Whittaker functions^ 


vr- _ C!\^2 ^ 

^ 1,2 ~ /- 

VC 

where 01,2 = T 


2r 

C(w) 

|A| 


A|2-w2 


(48) 

(49) 


Matching functions fi at r = R and using asymptotic 
Eqs. 65 j obtain a transcendental equation for the 
eigenenergies w of the sub-gap excited states: 


I 


V CXIai ,fi 

( ML,,, 


ai,ii 



FIG. 2: (Color online) Normalized energy of the first ex¬ 
cited state (^ = 1 and n = 0) as a function of the normal¬ 
ized island radius R/^ for different barrier transparencies A, 
BR^ — R^/lb <C 1, and according to Eq. (1161) . A = 20A 
corresponds to Ati ~ A, Ati ~ 0.75A, if A = 2A; and 
Ati « 0.54A, if A = A. 



FIG. 3: (Color online) Normalized energy split E+, Eq. (1371) . 
between edge and vortex Majorana fermions as a function of 
the normalized island radius R// for different barrier trans¬ 
parencies A = A and 2A, shift of the Fermi level U -I- SU — 
O.IA, and different magnetic fields '/BR — R/lb = 0.003 
(black) solid line, R/lb = 0.01 (red) dashed line, and R/h = 
0.03 (olive) dash-dotted line. 


In this expression the Whittaker functions Ma^^i{z) are 
taken at 2 ; = 2R/({uj), the Bessel functions Ja{z) at 
2 ; = ujR/v, and the prime means differentiation over z: 

= dMr, i,{z)/dz. An equation similar to E g. (f50l) 
was derived in Ref. and later corrected in Ref. [^. 

Each solution of Eq. (l50l) for the excited states can 
be characterized by a pair of quantum numbers: orbital 
number /r and principal number n. A similar classifica¬ 
tion scheme was used in Refs. lMlI^ . Our numerical anal¬ 
ysis shows that the lowest excited state, localized near the 
edge, corresponds to /i = 1 and n = 0. The energy of 





























the first excited state wqi as a function of R is plotted 
in Fig. [2] for different values of the barrier transparencies 
A/A. 

As one can see from this figure, the function a;oi(i?) de¬ 
creases when the island radius R increases. This function 
has a plateau at low R when the excited state approaches 
the edge of the continuous spectrum, that is, when wqi 
is close to Ati. The numerical results at large R/^ and 
A ^ A can be approximated by the ‘particle-in-the-box’ 
formula 


It follows from the results shown in Fig.j^that the energy 
gap between the first excited state and the Majorana 
fermion is of the order of the energy gap A in the SC if 
A > A and R < 10^. Thus, it is not necessary to have 
an ideal barrier between the TI and SC for a reasonable 
robustness of the Majorana fermion. For example, if the 
island radius is i? = 7.^ and A = 2A, we obtain wio ~ 
0.3A, which is much larger than the CdGM level spacing 
~ A'^/Ep for the vortex core states^^ (here Ep is the 
Fermi level of the SC). The SC island radius i? ^ 5^ is 
optimal for the stability of the Majorana state if A > A. 

In Fig. [3] we show the dependence of the energy 
split Ep, Eq. dSil), between edge and vortex Majorana 
fermions as a function of the island radius R/^ for dif¬ 
ferent barrier transparencies and magnetic fields. Com¬ 
parison of the results shown in Figs. [3] and [3] demon¬ 
strates that the splitting between two Majorana states is 
small {Ep A) and the Majorana state is rather stable 
(woi ~ 0.3 — 0.4A) under realistic values of parameters. 


V. DISCUSSION 

Below we discuss possible generalizations of our con¬ 
clusions beyond the constrains assumed in the previous 
sections, as well as connections of our results to that of 
other workers. 

We study the Majorana fermion near the edge of the 
s-superconductor island on the top of the topological in¬ 
sulator. In our case, the edge is a boundary between 
the 2D proximity-induced superconductivity and gapless 
surface of the topological insulator. 

One of the important parameters of the system studied 
is the island size. The applied magnetic field induces a 
vortex in the SC island and localizes a Majorana state 
near the island edge. However, if the island radius is 
comparable to the length scale ^(0) = v/X or ^ 5 c, the 
stability of the Majorana state deteriorates: first, due to 
the interaction of the edge and core Majorana fermions 
and, second, due to tunneling of the CdGM excited states 
to the edge. Therefore, in a small island there arises 
a rather peculiar picture of the GdGM statesi^ Thus, 
the condition (1351) is necessary for the existence of well- 
defined Majorana fermions in the system. 


An island of large radius may affect the distribution of 
the magnetic field in its vicinity. For our calculations we 
assumed that the magnetic field is uniform. This is true if 
R is much smaller than the effective London penetration 
depth Aeff in the SG island: R <C Aieff- However, even 
if this condition is violated, our results survive, provided 
that the magnetic field is not too strong 

BR^ = R?/ll < 1, (52) 


that is, the magnetic flux through the area of the island 
is smaller than the flux quantum. 

The latter inequality is the condition ensuring the va¬ 
lidity of our results. The main non-perturbative effect of 
the magnetic field is the stabilization of a vortex in the 
island, while the inhomogeneity of the magnetic field in 
the range r 4 may be studied perturbatively. Indeed, 
the generation of Majorana states at the vortex core and 
the edge depends on the parity of the vorticity quantum 
I, and is completely unaffected by the details of the mag¬ 
netic field distribution. 

Moreover, it can be shown that the relative correc¬ 
tion to the energy of the first excited state due to total 
screening of the magnetic field from the interior of the 
island is of the order oi R^/If. To evaluate such a cor¬ 
rection 5woi, we assume that the magnetic field vanishes 
for r < R. Then, following the procedure presented in 
the previous section we find 


SuJoi — 


R 

f 2TTr‘^dr ijj^iayijj 

y 0 _ 

2^2 oo 

^ f 2i:rdr 
0 



IfR 


^ Wqi, 


(53) 


if the conditions Eg. [52land A > A are valid. Eor a large 
island i? ^ this correction to the energy is small even 
if the magnetic field is quite strong R ^ Ij,. 

The case of stronger magnetic field or larger SG is¬ 
land, R/lb ^ 1, requires a separate consideration. How¬ 
ever, as before, Majorana fermion may exist only when 
the island hosts a vortex. This statement is quite natu¬ 
ral since the vorticity affects the quantization condition 
Eq. (I22p of fi. As a result, Majorana fermions can exist 
only if the number of vortices is odd^. In addition, for 
each Majorana fermion, its partner must exist because 
the fermion parity must be conserved. In the case of a 
singly-connected SC island the only possibility is to have 
one Majorana fermion near the island edge and another 
Majorana fermion near the vortex core, because two Ma¬ 
jorana fermions located in the same edge form a Dirac 
fermion. 

In Refs. 1^1^ the authors consider a semi-finite SC is¬ 
land on top of a TI in a transverse magnetic field. The 
Majorana fermion is delocalized at the edge between the 
SC and the TI. In such a geometry, the condition of 
single-valuedness of the wave functions Eq. (1331) becomes 
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a momentum quantization rule. If momentum quantiza¬ 
tion is ignored, then the Majorana fermion can exist at 
the edge of the SC even in the case of zero vorticity. This 
is just an unphysical artifact of using an infinite sample. 

One of the major driving forces behind the develop¬ 
ment of Majorana solid state research is the possibility of 
performing topological quantum computation. To be us¬ 
able in such a setup, the Majorana fermion must be well 
separated from non-topological excitations. We demon¬ 
strated that the energy of the first excited state localized 
at the edge could be as large as a fraction of the super¬ 
conducting gap, see Fig. [2] The time of an elementary 
braiding operation must be much less than the decoher¬ 
ence time^ caused by the hybridization of the edge and 
vortex Majorana fermions. In our system, such decoher¬ 
ence time increases exponentially with the radius of the 
island, see Eq. ( 1321 ). 

In addition to the excitations pinned at the island 
edge, there are non-Majorana states localized aX r > R. 
These are the Landau levels which we briefly discussed in 
Sec. IIIIBI Such states with large values of ^ are pushed 
away from the island by the centrifugal force. As a conse¬ 
quence, they cannot affect the Majorana state. At small 
/i, however, their wave functions can reach the island 
edge. Fortunately, the Landau levels are separated from 
each other by a gap of the order of Ab ~ v/Ib- Thus, 
at not-too-weak magnetic helds and finite U these states 
are shifted from the zero energy by an amount ^ Ab. 

Scanning tunneling microscopy (STM) might be a use¬ 
ful tool for investigating Majorana fermion o^^i^^ . In STM 
experiments, a Majorana fermion could manifest itself as 
a robust zero-bias peak. A detailed analysis of the STM 
spectroscopy of the edge Majorana fermions in the pres¬ 
ence of vortices is done in Ref. lHI. The stability of the 
zero-bias peak should be checked against variations of the 
chemical potential or the gate voltage, to distinguish the 
Majorana fermions from Dirac fermions, as discussed in 
Section IIIII 

When the magnetic field is varied, the strength of the 
zero-bias conductance on the island edge should oscillate 
when the number of vortices changes. These oscillations 
could be an additional proof of the existence of Majo¬ 
rana fermions. Observations of several vortices and mul¬ 
tivortices in superconducting Pb nano-islands were re- 
portedi^i^^ Thus, measuring zero-bias peak oscillations 
as a function of the vorticity is a realizable experimen¬ 
tal task. The coordinate dependence of |'i/'P for the edge 
and core Majorana states are shown in Fig. [4l The edge 
Majorana fermion penetrates in the island at the dis¬ 
tance ^ and outside the island at the distance 4- With 
the decrease of the magnetic field (and the growth of the 
magnetic length) the peak value in for the edge Ma¬ 
jorana state decreases and if R —>■ 0 the edge Majorana 
fermion becomes delocalized. Since the density of states 
is proportional to such a behavior can be observed 
as a zero-peak in STM measurement. Note, that the 
STM measurement could also reveal the zero-peak split¬ 
ting due to overlapping of the wavefunctions of the edge 



to the local density of states, as a function of radial position: 
(red) dashed line for the vortex Majorana fermion and (blue) 
dashed-dotted line for the edge Majorana fermion. The local 
density of states may be measured by the STM. (In the figure 
the value of for the edge Majorana fermion is multiplied 
by 100.) Solid line (green) line shows the proximity-induced 
superconducting order parameter A(r)/Ao, which vanishes 
when r > R, here Aq is the bulk value of A and R = h = 5^. 


and vortex Majorana fermions in small islands i? ~ ^ 
or due to the close localization of two superconducting 
islands. 

For numerical estimates, let us consider Tc = 10 K and 
A = 1.76Tc ~ 2 meV for a BCS-type superconductor. 
Assuming that A = 2A Ri 4 meV, then, A = A and ^ = 
2^. If we take the radius of the island R = 7^, then, the 
energy of the hrst excited state becomes wqi — 0.3A ~ 
5K. To evaluate the possible radius of the SC island we 
should estimate the value of the coherence length ^ = 
u/A, which depends on the Fermi velocity on the surface 
of the TI. In Ref. it was reported that v = 5.0 x 
10^ cm/s for the surface of Bi 2 Se 3 in vacuum, then, £ Ri 
200 nm for the value of the gap chosen here. In Ref. 0 it 
was obtained that v = 10^ cm/s on the interface between 
Bi 2 Te 3 and a nanoribbon, and then, f r; 40 nm. If so, 
then the appropriate value of R is of the order of several 
hundred nm. However, it has been reported in Ref. |4^ 
that the Fermi velocity on the surface of Bi 2 Te 3 can have 
a much lower value, v = 3 x 10® cm/s, and then, ^ r; 
1.2 nm and R could be of the order of several nm, which 
might be of the order of or lower than the coherence 
length ^sc in the bulk of the SC island. The case R < 
^SC is not optimal for the stability of the edge Majorana 
fermion because CdCM states in the vortex core in the 
bulk of the SC can affect the edge states. 

To conclude, we studied the electronic properties of 
a superconducting island in a magnetic field placed on 
the surface of a topological insulator. Majorana states 
arise only if a vortex with odd vorticity exists in the su¬ 
perconducting island. Non-topological excitations in our 
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structure are separated from the Majorana fermion by a 
significant gap, provided that the parameters are suitably 
chosen. A Majorana state may be detected in an STM 
experiment as a zero-bias peak, which is stable against 
variations of the gate voltage. The zero-bias conductance 
should oscillate as a function of the magnetic field. We 
here estimate the optimal parameters for the experimen¬ 
tal study of Majorana fermions in our system. 
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